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Abstract 

Using the SCRPA, we study the intersite interaction effect on the dynamics of N electrons system. 
We have considered an extended Hubbard model including intrasite and intersite interactions, and 
we have applied this model to a system of two neighbouring atoms containing a free electron. The 
application of SCRPA to this model allows us to study the intersite interaction effect on the ground 
state and the excitation energies of system. We show that the repulsive interaction between the 
electrons of the neighbouring atoms is the origin of an supplementary conductivity of the system. 
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1 Introduction 



The discovery of High- Temperature superconductivity has motivated a considerable effort in physics of 
strongly correlated electronic systems, and many theoretical models have been proposed[l]. The Hub- 
bard model [2, 3] is one of the simplest and more general description of an interacting electrons system 
on a lattice. In its simplest form, it includ the usual kinetic energy of electrons and the competition 
between the on-site electron-electron interaction. The standard Hubbard model is defined by the second 
quantized Hamiltonian [2, 3]: 

The first term of the eqs(l) represents the kinetic energy of electrons, and each electron has a possibility 
of hopping between different lattice sites. Cj t<T is the annihilation operator of the electron at a lattice site 
j with spin index a . c\ ^.is the creation operator of the electron at a lattice site i, so tij is the hopping 
integral from the site j to the site i. The second term represents the intrasite coloumb interaction with 
energy U, where ni j(T is the number operator of electrons at the site i with spin a . 

Recently, the Random Phase Approximation (RPA)[4], was used to resolve the standard Hubbard 
model [5, 6]. The RPA is an approach which treats seriously the correlations of system, and attempt 
to minimise the system energy. The standard Hamiltonian Hubbard must be developed as function 
of the creation and annihilation operators of the pair particle- hole(p-h), because our RPA regroups 
the electrons system on pair: particle-hole. The application of RPA to standard Hubbard gives non 
linear coupled equations, where the resolution is done by a SCRPA[7, 8]. In this paper, we consider an 
extended Hubbard model[9], where the intersite coulomb interactions are introduced. This model was 
shown to describe many interesting properties of high TC superconductors materials [10]. In this work, 
we apply the SCRPA method to this extended Hubbard model, and study the intersite interaction effect 
on the ground state and excitation energy of system. We show that the repulsive interaction between 
the electrons of the neighbouring atoms is the origin of a supplementary conductivity of the system. 



2 Extended Hubbard Model 

The standard Hubbard model with intrasite interaction explains some important physical phenomena 
like High- Temperature superconductivity [3], Mott-transition[ll]. To explain other physical phenom- 
ena observed in different areas of the solid state physics like magnetic and transport properties, it is 
convenient to take into account also the intersite interaction resulting from original coloumb repulsion 
modified by the polaronic effect. The extended Hubbard Hamiltonian is then given by [12, 13]: 

H = Yl *ij- c *>- c J> + r -X!"' • + \ Yl v ij~ ) - n h°- n J> + \ Yl v ij 2) - n i,v n j-° ( 2 ) 

U denotes the effective intrasite coloumb interaction, vjp ( v}P) describes the effective intersite coloumb 
interaction between the electrons in the lattice sites i and j, with the same spins (opposite spins). V^p 
and are not necessary equal. The model(2) cannot be solved in a general case. There is however, 

a special but non trivial case of finite number of sites, which possesses exact analytical solution[5]. 

In this work we will limit ourselves to a simple case, and will apply the general formalism of SCRPA to 
the two sites problem. We consider a closed chain in one dimension, with periodic boundary conditions 



N = 2. Our physical system is then equivalent to two neighbouring atoms containing a free electron. 
The Hamiltonian of the system is: 

H n = -t^Xc la .C2.,a + c 2a .c 1 ^) + U.{n 1A .n 1 . l +n2^.n 2 .x) (3) 

CT 

+V r i.y^ni |g .n 2 , g + V 2 .^2ni^.n 2 - a 

(7 cr 

where t = —t\i = —t^i- 

In order to apply the approximation SCRPA to the Hubbard model, it is necessary, first, to apply 
the Hartree-Fock approximation(HF) to the Hubbard model. In the HF method, we write the Hamil- 
tonian^) as function of quasiparticles operators, wich allow us to obtain the excitation spectrum of 
independent quasiparticles. The states \HF) are defined as: \HF) = a k . ^.a k . j_. \vac), where a k a is the 
annihilation operator of the mode (fc, cr); a k . a is related to c J;cr with the usual Fourier transformation: 

k is the momentum of state \k,a). The periodic boundary conditions suppose that cjy + j, a = c j,<y- 
With this condition, cqs(4) gives exp(— i. k .Rj) = 1, which have two solutions in the first Brillouin 
zone: fci = and k 2 = —n. Thus the Hamiltonian is then written as: 

H HF = E HF + {e-i.m^ + £2-n 2 ,cr} (5) 

This expression shows that in the HF approximation, the physical system has two possible states \HF) 
and \HF)* . \HF) (\HF)* ) is the Hartree-Fock ground state (excited state ), with the momentum: 
k\ = 0; below (and k 2 = —it; above ) the Fermi momentum, and the excitation energy: e\ (and 
e 2 ), where \HF) = a ki yd kl i • \vac) and \HF)* = ajj~ 2 r-a> k2 \vac). As in rcf [13], we define the HF 
quasiparticle operators by: bi iCr — a kl _ a and b 2M = a k2iCr ■ We have then b k ,a- \HF) = 0, for all k. 
In normal ordering of &i jCr and 62,0-; the Hamiltonian(3) becomes: 



where 



With 



H = H H f + -fffe=o + F[ k= -Tr (6) 



U + V 2 Vi -c-^ 2 

H k=0 = — - — (n k2A - n fcliT ) {n k2 ,i - n kui ) + — . ^ (h k2 , a - n kua ) 



H k=0 and H k =_n take into account the correlation between the number operators of the type: n ki . a n kj ^ a i 
in the ground state: k\ = (below the Fermi momentum ) and between the magnetic momentum 
operators of the type: f a -J~, in the excited state: k 2 = —it (above the Fermi momentum). 



3 Self Consistent Random Phase Approximation 



3.1 Formalism 

In order to apply the Formalism of SCRPA to the Hubbard model, it is convenient to use the particle- 
hole(ph-RPA) approximation, which regroup the physical system on pair. We can then define the RPA 
excitation operator as: 



Ql = - Vlh-bh.b p ) (7) 

p,h 

Where h (and p) are the momentum below (and above) the Fermi momentum. Eqs(7) shows that the 
excitation in the ph-RPA is done only by the creation or (annihilation) of pair: particle-hole via the 
operator bp.b 1 h (bh-b p ) with the amplitude x v ph (Up h )- The corresponding excited state of this excitation 
operator is \v) — Q* . \RPA), and the corresponding excitation energy is: 



E„ = 



(RPA\ 


Qv i 


H, Ql 


\RPA) 


(RPA\ [Q v ,Ql] | 


RPA) 



(8) 



Where \RPA) is the vacuum of this RPA excitation operator: Q v \RPA) 
The minimization of E v leads to usual RPA equations of type: 




v I V v I 

With the relations of the orthonormality conditions of the set |<3„; <5l|, we can express the elements 
of A and B by the RPA-amplitudes, and therefore we have a completely closed system of equations for 
amplitudes x and y. For our problem, we consider only the excitation operators, which conserve the 
spin, where the excitation is done only by the creation or annihilation of the pair: particle-hole with the 
same spin. 



Ql = x* v k\+xl.kl-ylkf -vlk; (9) 

With Kl = b\ a .b\ j^J\ - (M a ), K~ = b ha .b 2 ^/y/l - (M a ) and M a = n hrT +n 2 ,a, where the 
mean values (...) are taken with respect to the RPA vacuum (Q v \RPA) = 0). The SCRPA equation 
can then be written in the form: 











B 1l 








B n 


B ll 








-A TT 
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-A u J 



( ^ \ 

-I 

V vl J 



Where the SCRPA matrix elements are given by: 



( X T \ 

V vl J 



(10) 



and 



B aiTl = {[K-,[H,K-]]) 
The orthonormality relations of the set ^Q v ; Q f v |, give: 



^TT = A U =A A n =A l1= A' 

S TT = B ll = B B n = Bit = B' 
From the Hamiltonian given in eqs(10), we can writte the SCRPA matrix elements as: A = B + 21, 
A 1 = B' , where: 



A = 2.t+(U 



Vi 



l 



2 I 1 - (M a ) 



(x a .x a + y a .y a + 2.x a .y a ) 



and 



where 



[/- V 2 1 



(Ma) 



2 1 - (M a ) 

2 E l#l 2 

1 + 2 E \VZ\ 2 



(11) 
(12) 

(13) 



so, the ph-RPA matrix takes the form: 



A - 21 A' \ 



(A A 1 

A 1 A A 1 A-2.t 

2.t-A -A' -A -A' 

\ -A' 21 -A -A' -A ) 

This ph-RPA matrix has two positive roots: 



( x\\ 

V vl J 



V vl J 



si = 21. 



A- A' 
t 



1 and e 2 = 21. 



A + A' 
t 



1 



(14) 



The corresponding eigenvectors are: V\ = [xi.,—Xi,yi,—yi] and V 2 — [x 2 ., — x 2 , y 2 , — y 2 ], respectively. 
Where 



Xi 



A-A' + ei 
' A-A'-2.t 



X 2 



A + A' + ei 
' A + A' -21 



■yi 



and 



yi = -\/2- 



A- A' + ei 
A-A>- 21 



- 2 



V2 = -\2. 



A + A'+e 2 



- 2 



(15) 



(16) 



A + A' - 21 / 

So, like the HF approximation, in ph-RPA, our system have tow excitation energies eiand s 2 , but they 
are coupled. Thus in this work, we solve a system of a coupled equations numerically by iteration leading 
to a SCRPQ solution which are quasi identical to the exact result. 



Figure 1: 



3.2 Results and discussion 

To show the effect of intersite interaction on the energy of the system, we have studied th evolution of 
the ground state and excited energies in term of the interaction V\ and V2. 

In figurel we plot the variation of the ground state energy Escrpa = (0 \H\ 0) as function of the 
two parameters of the intersite interaction V\ and V2 ■ 

The result shows that the ground state energy is quasi independent on V± but, the variation of 
Escrpa become more important when we introduce the intersite interaction with the opposite spins: 
V2. This results can be explained by the fact that the SCRPA include only, for the fundamental state, 
the correlations between the particles with different spins: |0) = ( cj + c\ jt . J \HF) 



Thus the only type of interaction wich is of interest is V2 ■ In the following we analyse the V2 effect 
on the dynamics of system. In figure 2 and 3 we plot the variation of the gorund state energy and 
the excitation energies of SCRPA, respectively, as function of the repulsive intrasite interaction U for 
different values of the intersite interaction V2. 




Figure 2: 



Figure 3: 



The results show that, when U increases, the gap A = e 2 — t\ between the two excitation energies 
increases too, and so, the jump of electrons between these two states becomes more difficult. Thus, we 
can conclude, that for a fixed value of U, while the intersite interaction is repulsive (attractive), the gap 
A becomes less (more) important. These remarks allow us to assume that repulsive interaction between 
the electrons of the neighbouring atoms is the origin of supplementary conductivity of the system. 

4 Conclusion 

In this paper, the SCRPA approximation was used to solve the extended Hubbard model given in 
Eqs(3). The quality of the SCRPA method has been investigated in a previous work by Jemai[5], in 
which he has shown a remarkable agreement between SCRPA method and excact results for the standard 
Hubbard model. In our work, we have extended this technic to study the intersite interaction effects 
on the dynamics of the electrons in the two sites with (n^f) = (n^i). We have shown that the gap 
between the excitation energies: (A = e 2 — £i) are correlated with the intersite interaction energy Vi. 
This result allow us to suppose that the repulsive intersite interaction (between the electrons of the 
neighbouring atoms) is the origin of a supplementary conductivity of the system. In future work[14], we 
propose to solve the 4-sites case (plaquette), which may be very important for the explanation of high 
T c superconductivity, by considering the many plaquette configurations in 2D. 
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